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Abstract 

The Wigner function shares several properties with classical dis- 
tribution functions on phase space, but is not positive-definite. The 
integral of the Wigner function over a given region of phase space can 
therefore lie outside the interval [0, 1] . The problem of finding best- 
possible upper and lower bounds for a given region is the problem 
of finding the greatest and least eigenvalues of an associated Hermi- 
tian operator. Exactly solvable examples are described, and possible 
extensions are indicated. 



1 Introduction 

Since its introduction £Q in 1932, the Wigner distribution function has proved 
an important tool in many areas of quantum physics and chemistry. Recently, 
in quantum tomography [2] , the Wigner function has been reconstructed from 
experimentally measured probabilities associated with a quantum system 
in a definite state. Given a finite amount of data, it is clear that such a 
reconstruction must inevitably be limited to a finite region of phase space. 
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And here there arises a problem peculiar to quantum tomography. In the 
classical case, if the probability density has been reconstructed on a finite 
region, and if its integral over that region, the 'probability mass' on the 
region, is close to 1, then it is certain that the density on the remainder of 
phase space is small, in the sense that the probability mass there must be 
small. But in the quantum case, because the Wigner function takes negative 
as well as positive values in general, the 'quasiprobability mass' on regions of 
phase space outside one particular region may be large positive or negative, 
even when the quasiprobability mass on that particular region is known to 
be close to 1. This is the primary motivation for our study, which aims to 
provide best possible upper and lower bounds on the integral of the Wigner 
function over any given region S of phase space. 

In previous work it has been shown that this is equivalent to finding 
the greatest and least eigenvalues of a corresponding Hermitian operator S, 
which we refer to as the 'region operator' for S. It is an operator analogue 
of the classical characteristic function 5$, which has the value 1 for points 
inside S, and for points outside. Here we describe some exact results and 
the basis on which they can be obtained, and indicate some directions of 
possible generalization. 

2 Mathematical background and results 

For a system with one degree of freedom, the invertible mapping from func- 
tions A(q,p) on the phase plane T to linear operators A on Hilbert space is 
accomplished with the help of Weyl's kernel operator 



Here Tr denotes the trace, and q, p are canonical operators. The generaliza- 
tions to many degrees of freedom are obvious. From this point we work with 
dimensionless variables, in effect setting h = 1. 

As special cases, we have the Wigner function W p (q,p) coresponding to 
the operator p/(2n), where p is the density operator defining the state of the 
quantum system; and the region operator S, corresponding to the character- 
istic function Ss(q,p) of a given 2-dimensional region S G T : 




W p (g,p) = Tr(W(q,p)p), S 



<*s(q,p)W(q,p) dqdp 



L 



W(q, p) dqdp . 

(2) 
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In the last step we have used the defining property of the characteristic 
function. 

Given a smooth contour C G T, it is natural to define also an associated 
'contour operator' C by 

6 = Jw(q,p)S c (q,p)dqdp. (3) 

Here 5c(q,p) is the 'characteristic (generalised) function' for the contour, 
defined by the property that, for every smooth function FonT, 

J F(q, p)S c {q, p)dqdp = J F dl , (4) 

where dl the element of length on C . Then we can write 

6=f Wdl. (5) 
Jc 

For example, if the contour is defined by the graph of a simple function, such 
as p = F(q), a < q < b, then Sc(q,p) = 5{p — F(q)) for a < q < b, and 
$c(q,p) — elsewhere. Going one dimension further down, we can define the 
'point operator' 

P = J 6 P (q,p)W{q,p) dqdp = W(q ,p ) . (6) 

corresponding to the point (qo,Po) G T. The characteristic function for the 
point, Sp(q,p) = S(q — qo)8(p — Po), is more singular than that for a contour 
C, which is more singular than that for a 2-dimensional region S in the 
plane. Clearly these notions can be generalized to regions R of dimensions 
2f, If — 1, . . ., in the 2/-dimensional phase space of a system with / degrees 
of freedom, with the introduction of corresponding characteristic functions 
5r of increasing singularity, and corresponding region operators R. 

Critical to what follows is Wigner's observation pP that, when the quan- 
tum system is in the state p, the mean value of any operator (such as R) is 
given in terms of the corresponding phase space function (in this case 5r) 
and the Wigner function, by 

(R)= [ 5 R W p dT= I W p dV, (7) 
Jr Jr 

where dV is the volume element in R. Since the expectation value (R) must 
lie between the greatest and least eigenvalues of R (or more generally, between 
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the best possible upper and lower bounds on its spectrum), the same must be 
true of the last integral in (JJJ). Moreover, the bounds are attained (or more 
generally, approached arbitrarily closely) if p is a pure state corresponding 
to the associated (generalized) eigenvector of R. In short, if we can solve the 
eigenvalue problem for R, we have best possible bounds on the last integral 
in (|7ji. and we know when and how those bounds are attained. 

The circular disk of radius a, centred on the origin in the phase plane 
is invariant under canonical transformations generated by the classical har- 
monic oscillator Hamiltonian H = p 2 + q 2 , and it follows that the corre- 
sponding region operator S a commutes with the quantum oscillator Hamil- 
tonian operator, and has the same eigenvectors, namely the number states 
|n), n = 0, 1, .... This is significant for quantum tomography, where such 
states are amongst the easiest to measure 0. The eigenvalue of S a on \n) is 

r 2 

X n (a) = (-l) n L n (2u)e- U du } (8) 
Jo 

where L n is the Laguerre polynomial. For a given a, it is then possible 
to see which eigenvalue is smallest (say Aat(o)), and which is largest (always 
A (a)), and hence what are the lower and upper bounds on the integral of the 
Wigner function over the disk. These bounds are attained when the system 
is in the pure state \N) or |0), respectively. The eigenvalues of the contour 
operator C a corresponding to the boundary of the circular disk are, again on 
the number states, /i«(a) = d X n (a)/da. There is a remarkable Lie- algebraic 
structure [3] associated with these spectra. 

For the contour operator Cl corresponding to a straight line segment of 
length L in the phase plane, say the segment along the p-axis from p = —L/2 
to p = L/2, we find 6\ a continuous spectrum and generalized eigenfunctions: 

„(*) = ±^ MkL) , (0 < k < co) ; = -j= { . (9) 

This segment is invariant under canonical transformations generated by the 
free-particle Hamiltonian H = p 2 , and also under the parity transformation 
which sends q, p to —q, —p, hence the appearance of the common eigenfunc- 
tions of the quantum free-particle Hamiltonian operator and parity operator 
in ©. 

From these simple results for disks, circles and straight segments we can 
build up results for more complicated regions in phase spaces of higher di- 
mensions. For example, in the 4-dimensional phase space of a 2-dimensional 
quantum system, we can find best possible upper and lower bounds on 
the integral of the Wigner function on the rectangle —L\/2 < q\ < L\/2, 
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— L 2 < q 2 < L 2 /2 within the 'Poincare section' p\ = p 2 = 0; or on the 
cylinder qf+pf < a 2 , — L/2 < p 2 < L/2, q 2 = 0; or on its surface, and so on. 

Investigations continue towards a "spectral theory" of quasiprobability 
mass for bosonic systems, and into analogous results for spin and other fi- 
nite dimensional quantum systems, with expected applications in quantum 
tomography. 

Notes 
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